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Abstract: HF radar systems are widely and routinely used for the measurement of ocean 
surface currents and waves. Analysis methods presently in use are based on the assumption 
of infinite water depth, and may therefore be inadequate close to shore where the radar 
echo is strongest. In this paper, we treat the situation when the radar echo is returned from 
ocean waves that interact with the ocean floor. Simulations are described which 
demonstrate the effect of shallow water on radar sea-echo. These are used to investigate 
limits on the existing theory and to define water depths at which shallow-water effects 
become significant. The second-order spectral energy increases relative to the first-order as 
the water depth decreases, resulting in spectral saturation when the waveheight exceeds a 
limit defined by the radar transmit frequency. This effect is particularly marked for lower 
radar transmit frequencies. The saturation limit on waveheight is less for shallow water. 
Shallow water affects second-order spectra (which gives wave information) far more than 
first-order (which gives information on current velocities), the latter being significantly 
affected only for the lowest radar transmit frequencies for extremely shallow water. We 
describe analysis of radar echo from shallow water measured by a Rutgers University HF 
radar system to give ocean wave spectral estimates. Radar-derived wave height, period and 
direction are compared with simultaneous shallow- water in- situ measurements. 
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1. Introduction 

HF radar systems are widely used internationally to provide continuous monitoring of ocean waves 
and currents for a large range of environmental conditions. 

Within the US, coastal ocean current mapping with HF radar has matured to the point where it is 
now considered an important component of regional ocean observing systems. A mid- Atlantic HF 
radar network now provides high resolution coverage within five localized networks, which are linked 
together to cover the full range of the mid-Atlantic coastal ecosystem. Similar regional networks 
around the US coastline are being formed into a national HF radar network. 

While much of the focus of these networks until now has been on offshore current mapping 
observations, a longer-term objective is to develop and evaluate near-shore measures of waves and 
currents. These investigations aim to understand the interaction of waves in the shallow coastal waters 
and how energy is transformed into the creation of dangerous rip currents along the New-Jersey/Long- 
Island shorelines. Rutgers University radars cover these coastal regions at multiple frequencies from 
4.5 to 25 MHz. Their echoes contain information on both currents and waves from deep water up into 
the shallow coastal zone, providing an excellent archive for such studies. This paper describes the 
analysis of both simulated and measured radar echo to demonstrate the effect of shallow water on radar 
observations and their interpretation. 

Radar sea-echo spectra consist of dominant first-order peaks surrounded with lower-energy second- 
order structure. Analysis methods presently in use assume that the waves do not interact with the 
ocean floor, see [1, 2, 3] for phased-array-antenna beam-forming systems; and [4] for systems with 
compact crossed-loop direction-finding antennas, such as the SeaSonde. 

The assumption of deep water is often invalid close to the coast and for broad continental shelves, 
and is particularly inadequate to describe the second-order sea-echo used to give information on ocean 
waves., as second-order echo is often visible above the noise only for close ranges. To interpret this 
echo correctly, we show that the effects of shallow water must be taken into consideration. 

In Section 2, we give the basic equations describing radar echo from shallow water, expanding on 
the previous description given in [5]. In Section 3, simulations are used to illustrate the effects of 
shallow water on waveheight, Doppler shifts and spectral amplitudes in radar sea-echo spectra, to 
investigate limits on the existing theory and to define depth limits at which shallow-water effects must 
be included in the analysis. The effects of shallow water on the radar spectrum are illustrated using 
measured spectra. In Section 4, methods are applied to the interpretation of measured radar echo from a 
Rutgers University radar to produce wave directional spectral estimates, which are compared with 
wave observations from a bottom-mounted Acoustic Doppler Current Profiler (ADCP) moored in the 
second radar range cell. 

2. Radar spectral theory 

It follows from the solution of the equations of motion and continuity that long ocean waves are more 
affected by shallow water. We define the depth at which waves interact with the ocean floor by the 
approximate relation: 
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d/L < 1/8 (1) 

where d is the water depth and L is the dominant ocean wavelength. The deep-water analysis must 
be modified to allow for shallow-water effects in the coupling coefficients, the dispersion equation 
refractive effects on wave direction, and the directional ocean wave spectrum itself. We only consider 
water of sufficient depth that effects of wave energy dissipation such as breaking and bottom friction 
may be ignored; thus we operate in the linear wave transformation regime. As a general rule, this 
assumption is valid when the water depth is greater than 5% of the deep-water wavelength. 

Applying the lowest-order shallow-water dispersion equation to first-order backscatter from the sea 
gives the following equations for k] , the first-order spatial wave vector and co\ , the temporal 

wavenumber of the ocean waves in shallow water producing the backscatter. In this document, a 
subscript or superscript s indicates a shallow-water variable; its absence indicates a deep-water 
variable. 

K s — - k q 



o) 1 = m x m n 

s B 



(2) 



where k 0 is the radar wave vector, of magnitude k 0 , and C ° B is the Bragg resonant frequency in 
shallow water which is given by: 

co B - <s]2gk 0 tanh(2£ 0 <i) (3) 
with g the gravitational constant. The analogous relations for second-order backscatter are: 

co 2 = mJgk s tarih(k s d) - m 1 J^k x<\nh(k s d) 

(4) 

k k 

where v? v are the spatial wavevectors (with magnitudes k s , k s ) of the two shallow- water, first- 
order ocean waves interacting to produce the second-order backscatter. m, m'are equal to +1, -1 for 
waves moving toward, away from the radar respectively. 

The electromagnetic coupling coefficient has the same form as for deep water [5] but with shallow- 
water wavevectors: 



I'm; =0.5 



{k,.k 0 )(k s .k 0 )ikl-2k s .k s 



(5) 
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where A is the normalized surface impedance. The hydrodynamic coupling coefficient, derived by 
Barrick and Lipa [6] through solution of the equations of motion and continuity, is a function of water 
depth: 



P 



k + k 



[kk'-k s .k s ) 



nun \ kk 



or + co 



CO" - COl 



CO 



j( mjffof csclr (* A d) + (///' JxiT ) 3 csclrU'c/) 



(6) 



where k and k' are the spatial wavenumbers of the scattering waves in deep water. The deep- and 
shallow- water spatial wavenumbers are related as follows: 



k = k tmh(kd) k } = k\ tanh(^ ' d) 



(7) 



The total radar coupling coefficient r s i s the coherent sum of the hydrodynamic and 
electromagnetic terms 



1 1 EM T L H 



(8) 



It can be shown from these equations that at constant wavenumber, the coupling coefficient 
increases as the water depth decreases, resulting in an increasing ratio of second- to first-order energy 
as the depth decreases. 

In the following analysis, we assume that the deep-water directional wave spectrum is spatially 
homogeneous and that any inhomogeneity in shallow water arises from wave refraction. When energy 
dissipation can be neglected, it follows from linear wave theory that since the total energy of the 
wavefield, is conserved, the shallow-water wave spectrum expressed in the appropriate variables is 
equal to the deep-water spectrum [7] : 



S s (k s ) = S(k ) 



(9) 



where the deep- and shallow-water wave vectors are related by Snell's law and the dispersion 
equation: 



k cos(0 + J3) = k s cos(0 s + JS) 



(10) 



k = k s tanh(k s d) 



(11) 
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Here r is the angle between the radar beam and the depth contour and * ' are the angles between 
the radar beam and the shallow-, deep-water ocean waves respectively. Figure 1 illustrates refraction at 
a contour between regions of differing depth. 

Figure 1. Schematic geometry of the radar beam and an ocean wave train at a depth 
contour, denoted by the dashed line. Wave angles are measured counter-clockwise from 
the radar beam to the direction the wave is moving. Increasing * ' by 180° would 
define an incoming wave. 




The shallow- and deep-water rms waveheights are given by: 
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oo oo 



Substituting (10) and (11) into (12) gives the following relations which are useful for deriving the 
shallow- from the deep-water wave spectrum and vice versa: 



Rl = \ \ S s(KA)J(KA)kdk s d0 s H] = \\ S{k ,0)J-\k s9 0 s )k s dk d0 (13) 

0 0 

where the Jacobian ^^K^s) is given by: 



o o 



fdk) 
Id* J 



(do) 

fdk) (do) 
[doj 0 [del, 



s y 9 



fdk) fde) 
{dkJXdoJ, 



S y 0 x S' k 



1 + 



k s d sec h 2 (k s d) 
tarih(kd) 



sin(^+^) 
sin(0 + ]B) 



(14) 



The first- and second-order radar cross sections in shallow water at frequency CO and azimuth 
angle q> are given by: 



a\{a),(p)=kl S s (2k 0 ,<p + (m'+l)-)S(a)-m'a) B ) 



(15) 



n'=±l 



where S s (k, a) ) is the directional ocean wave spectrum for wavenumber k and direction a □ . 

2 TV oo 

cT*(co,<p)=k 4 0 X $ j\r 2 s \S s (k s ,0 s + <p+mx) 

in, in '=±1 o -oo 

• S s (k' s ,0 s + (p+m'ri)S(a)- m^Jgk s t<mh(k s d) -m^gk s t<mh(k s d) ^ s dk s dO s 



(16) 



where the coupling coefficient T s is given by (8). The values of m and m' in (16) define the four 
possible combinations of direction of the two scattering waves. Common numerical multiplicative 
constants in (15) and (16) have been omitted. It can be shown from (4) that the wavenumbers of the 
scattering waves are related as follows: 



k\=^k* + 2k s cos(0 s ) + l 



(17) 



To compute the second-order integral in (16), we choose as integration variables s and the deep- 
water angle 0 . In terms of these variables (16) becomes 
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m,m'=±l o 



3^ 



dhdO 



(18) 



where 



and 



, #) = m^Jgk -m'yjgk s tanh(k s d) 



l(k s ,0) = F 2 S S(k,0+<p+mx)S(k',0+<p+m'x)k 



and where we have substituted (9) for the shallow water directional spectra. The factors 
— — are obtained by differentiation using (10), (1 1) and (19). 

Ok) n 



(19) 



(20) 



k and 



To calculate the integral in (18), it is first reduced to a single-dimensioned integral using the delta 
function constraint. The remaining integral is computed numerically. 
Frequency contours are defined by: 



ct)-h(k,0) = O 



(21) 



which is solved for k as a function of ^for a given value of 0) . Due to wave refraction, the shallow 
water angle and wavenumber have discontinuities when the deep-water wave moves parallel to the 
depth contour, i.e. when 



0 = -j3, Tt-P 



(22) 



where /3 is the angle between the radar beam and the depth contour. 

Frequency contours are hence also discontinuous due to this effect at deep-water wave angles 
defined by (22). Examples of frequency contours for deep- and shallow-water are shown in Figure 2, 
plotted in normalized deep-water spatial wavevector space k/(2k 0 ). Normalized components p, q are 

defined so that p is along the radar beam and q perpendicular: 



p = (k o + kcos(0))/(2k o ) 
q = ksin(0)/(2k o ) 



(23) 



The discontinuities in the frequency contours are more pronounced when the contour is drawn in 
shallow-water wavenumber space, as it follows from (10), (11) that there are discontinuities in the 
shallow-water wave angle due to wave refraction. 
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Figure 2. Examples of frequency contours for water of depth 10m (continuous lines) 
compared with the corresponding contours for deep water (dashed lines). 
Normalized frequency: 0)Io) b=\2, >^=60deg. 
Radar frequency: (a) 5Mhz , (b) 25Mhz 
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It can be seen from Figure 2 that the deep-water ocean wave numbers corresponding to a given radar 
spectral frequency change with depth: they become either greater or smaller than the deep-water 
values, depending on the wave direction. This results in the frequency of second-order peaks in the 
radar spectrum changing with water depth. 

The effects of shallow- water on measured radar spectra are illustrated in Figure 3, which shows 
measured spectra from a 5 MHz radar in five radar range cells, with distances ranging from 18km to 
60km. As the water depth decreases, the second-order energy increases relative to the first-order and 
the frequency displacement between the first- and second-order peaks decreases. In the outer ranges, 
the second-order structure is almost the same from range cell to range cell, as the water is effectively 
infinitely deep. 
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Figure 3. Spectra from a 5MHz SeaSonde monopole antenna. Range/ Water depth: 
(a) 18km/ 5 -20m (b) 30km/10-50m (c) 42km/ 20-70m (d) 48km/ 35-80m (e) 54km/ 
40-100m 
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3. Narrow-beam radar spectral simulations 



To gain insight into the effects of shallow water, simulated radar echo spectra were calculated for a 
narrow-beam radar, using the model directional wave spectrum defined in [8] which consists of the 



Sensors 2008, 8 



4620 



sum of two terms: a continuous high-frequency wind wave spectrum and a swell component that is an 
impulse function in both wavenumber and direction. The swell component is defined by 



s(k,e s ) = H s 2 S(0-0 s )S(k-k s ) 



(24) 



where s , s ' * are the specified rms waveheight, direction and wavenumber. For this model, four 
sharp spikes occur in the radar spectrum. Here we consider only the second-order sideband for 
which m=l, m'=l , and assume water depths in the range 5- 100m and radar transmit frequencies of 
5Mhz and 25Mhz. For these values, it can be shown numerically that Doppler frequencies are always 
greater than the positive Bragg frequency. The radar beam is taken to be pointing perpendicular to 
parallel depth contours (i.e @ = 90° in Figure 1) 

3.1 Effect of water depth on waveheight 

For our model it follows from (13) that the relationship between the shallow- and deep-water rms 
waveheights is given by: 



H: = H m 



sin(^+/?)/sin(^+/?) 
I [ tanh(fc*<i) + k*d sec h 2 (k*d)~\ 



(25) 



This relationship is of course independent of radar frequency and has many angle symmetries. 
Figure 4 shows the ratio plotted as a function of depth for different wave directions. 

Figure 4. The ratio of shallow- to deep water waveheight plotted vs. depth for a 12 s 
wave. Wave direction in deep water relative to the radar beam: Red 180°, Blue 135° 
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It can be seen from Figure 4 that the waveheight initially decreases with decreasing depth as the 
wave enters shallow water but increases at depths below about 20m, which agrees with [7]. 



3.2 Effect of water depth on Doppler shifts 



It follows from (3) that for a given radar frequency, the Bragg frequency decreases with depth, 
causing the Bragg peaks to move slightly closer together. Figure 5 shows the Bragg frequency plotted 
as a function of depth. 



Figure 5. Bragg frequency plotted as a function of depth. 
Radar transmit frequency: Red 5Mhz, Blue 25Mhz 
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It can be seen from Figure 5 that the change in the Bragg frequency with depth is small. 

Figure 6 shows the displacement of the second-order peak from the Bragg frequency plotted as a 
function of depth for an 1 Is wave moving at different angles with respect to the radar beam. 

It can be seen from Figure 6 that as the water depth decreases, the second-order peak shifts toward 
the Bragg frequency for waves moving toward the radar, and further away for waves moving away 
from the radar. This is consistent with the two branches of the contour plot as shown in Figure 2. This 
effect is more marked for lower radar frequencies and can be seen in the measured spectra shown in 
Figure 3 in which the second-order peak moves closer to the first-order as the range from the radar and 
water depth decrease, with waves moving toward the radar. 
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Figure 6. The frequency shift of the second-order peak from the Bragg frequency for an 
lis wave, (a) 5Mhz (b) 25Mhz. Angle between wave and radar beam: Yellow 0°, Blue 
45°, Green 135°, Red 180° 



BP 0 2 
X 

■I 0.18 



« o.ie 

? 

§ 0.14 



1 0.12 



0 1 



ts 
c 
o 



c 0.08 
a? 



aoe 



0.04 



aj 0.02 

ST 

m 



























L 




(a) 












1 



J 




^ B fl j 

L j 


i f 


r i 


i j 


r 


r i 




r 

-— - 


r i 



= — 


r i 



j 


r 






















































■ 












tap _ J 


j 


j 

_*. 


■ i 




_____ j 


























______ 




______ 







9 10 11 

Depth (m) 



12 13 



14 15 



m 0.2 

x 

0.18 

I 

a 0.16 

? 

§ 0.14 



1 0.12 



0 1 

0.08 
006 
; 0.04 



y 

§ 0.02 
£ 0 































(b) 












1 




r 

1 - - - ^ 

L 




r 


r 



L J 


r l 



L 


r i 



L J 


r i 



L J 




r 



It 




















































1 " 


















f j 


1 


1 


1 




























______ 



















9 10 11 

Depth (m) 



12 13 



14 



15 



3.3 Effect of water depth on radar spectral amplitudes 

It is shown in [8] that for the impulse-function model defined by (24), the ratio R of the second- 
order to first-order energy is given by: 
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R = 2H* s 2 \F 2 s \ (26) 

ft fc* 

where the coupling coefficient T s is evaluated at wavevectors defined by * ' * . T 5 increases with 
decreasing depth and increasing wave period at a given radar frequency as illustrated in Figures 7 and 
8, which also show that shallow water has a greater effect as the radar transmit frequency decreases. 

Figure 7. The absolute value of the coupling coefficient T s vs. depth for a 9 sec wave. 
Radar frequency: Red: 5Mhz, Blue: 25Mhz. 
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Figure 8. The absolute value of the coupling coefficient vs. depth for waves of different 
period. Radar transmit frequency: 5Mhz. Wave period: Red 15s, Blue 12s. Green 9s. 
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Since the coupling coefficient increases as the depth decreases, it follows from (26) that the second- 
order energy will increase with respect to the to first-order. This effect can be seen in the measured 
radar spectra shown in Figure 3. Figure 9 shows the theoretical ratio of the second- to the first-order 
energy obtained from (26) using our model for an 1 Is wave. 

Figure 9. Ratio of second - to first-order energy for an lis wave. Significant 

waveheight: 2.4m. Radar transmit frequency: (a) 5 Mhz, (b) 25 Mhz. 

Angle between wave and radar beam: Yellow 0°, Blue 45°, Green 135°, Red 180° 
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It can be seen from Figure 9 that the ratio of the second- to the first-order energy exceeds unity (i.e. 
the calculated second-order energy exceeds the first-order energy) for depths less than about 8 m for a 
5MHz transmit frequency and for depths less than about 10m for a 25 MHz transmit frequency. 

This subsection demonstrates an important point. Since we have shown that the waveheight itself 
actually decreases slightly upon moving into shallow water, while the second-order echo increases 
significantly due to the rapid growth of the coupling coefficient, wrongly using deep-water inversion 
theory to estimate waveheight will overestimate this important quantity. We note that all previous 
treatments and demonstrations of wave extraction have been based on deep-water theory, even when in 
fact many of the radar observations have been made in shallow water. 

3.4 Effect of water depth on breakdown of theoretical model 

When the magnitude of the second-order energy approaches that of the first-order, it is apparent that 
the perturbation expansions on which (15) and (16) are based are failing to converge and they therefore 
cannot provide an adequate description of the radar echo. This effect is similar to the well known radar 
spectral saturation occurring when the waveheight exceeds a limit defined by the radar transmit 
frequency. Above this waveheight limit, the radar spectrum loses its definitive shape and the 
perturbation expansions fail to converge. The deep-water saturation limit on the significant waveheight 
W Sat (defined to be four times the rms waveheight) is given approximately by the relation: 

W Sat = 2/k (27) 

For shallow-water, the saturation of the radar spectrum is exacerbated by the increase of the 
coupling coefficient and the radar spectrum saturates for waveheights less than that defined by (27). 
We here define the shallow-water saturation limit W s s at for the model to be that waveheight for which 

the second-order energy equals the first-order, and the ratio R is given by: 

R = 1 (28) 

In practice the theory may fail before this limit is reached. W Sat and W$ at are plotted vs. depth in 
Figure 10 for two different radar frequencies. At depths of 30m the saturation limits are approximately 
equal. At depths less that 30m, the shallow-water limit drops off sharply, particularly for the lower 
transmit frequency. Thus the radar spectrum can be expected to saturate at lower values of waveheight 
in shallow water. 

For waveheights above the saturation limit, the waveheight predicted by the theory will be too high. 
However the theory cannot be applied at all when the second-order spectrum merges with the first, as 
then separation is not possible. 
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Figure 10. Significant waveheight saturation limits for an 11 -second wave coming 
straight down the radar beam. Radar transmit frequency: (a) 5 Mhz , (b) 25 Mhz 
Red: deep-water saturation limit W Sat , Blue: shallow-water saturation limit W$ at 





3.5 Depth limits for significant shallow-water effects 

We estimate depths for which shallow- water effects become significant as follows: For first-order 
echo, the depth limit is defined by equality in (1). At this depth, the Bragg frequency defined by (3) is 
96% of its deep-water value. For second-order echo, we define the depth limit D s at which shallow- 
water effects become significant as the value at which the coupling coefficient defined by (8) exceeds 
1.25 times the deep-water value. Figure 11 plots the depths D s vs radar transmit frequency for an lis 

wave. 
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Figure 11. Depths at which shallow-water effects become significant vs. radar transmit 
frequency. Red: second-order echo. Blue: First-order echo. 




Radar transmit frequency (MHz) 

Figures 10 and 11 help in assessing the validity of the existing deep-water methods. However 
they are based on a wave model (24), which is quite restrictive: waves of a single wavelength are 
assumed to come down the radar beam. Also Figure 11 applies only to an lis wave. Performing 
similar studies for more general wave spectral models is beyond the scope of this paper. However, we 
observe that: (a) Shallow water effects are stronger for longer ocean wavelengths (b) Second-order 
radar spectra for m=l, m'=l are strongest for waves down the radar beam, (c) The stronger the second- 
order energy for a given waveheight, the sooner the radar spectrum will saturate as waveheight 
increases. Therefore shallow- water effects will be more marked at a given waveheight for a broad 
nondirectional spectrum that includes longer wavelengths e.g. the Pierson Moskowitz model (32). 
These differences would probably not be large however, due to the sharp cutoff of wave-spectral 
models for long wavelengths. The opposite effects would be expected for spectra that include wave 
directions not directly down the radar beam e.g. a cardioid directional distribution. To summarize these 
effects: W s s at will be less and D s will be greater than the values shown in Figs. 10 and 11 for the 

following changes from the wave spectrum (24): broad nondirectional spectrum, wave period > lis. 
W$ at will be greater and D s will be less for broad directional distributions, waves nonparallel to radar 

beam , wave period < 1 Is. 
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4. Application to measured data 



4.1 Data set 



The results presented here are based on analysis of 10-minute radar spectra measured by a 25MHz 
SeaSonde located at Breezy Point, NJ. The time period from December 29 to 30, 2005, was chosen 
because simultaneous coverage provided by the SeaSonde and a bottom-mounted ADCP allowed a 
direct comparison to be made between results from the two sensors. The ADCP was located in the 
second radar range cell in water of depth 8m. The bathymetry in the area and the locations of the two 
sensors are shown in Figure 12. 

Figure 12. The coastline and bathymetry (contours in meters) around Breezy Point, 
New Jersey, showing the positions of the SeaSonde and the bottom-mounted ADCP. 




In our analysis, depth contours near the radar are assumed to be parallel to shore and the depth 
profile is obtained from Figure 12. 

Figure 13 shows measured spectra from the Breezy Point SeaSonde at three ranges: the second- 
order energy can be seen to increase relative to the first-order as the water depth decreases. 
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Figure 13. (a) Spectra measured by the 25MHz SeaSonde at Breezy Point, at 1 :00pm 
12/30/2005. Range: (a) 3 km (b) 6 km (c) 9 km. 
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4.2 Interpretation of the radar spectra 

Lipa and Barrick [5] describe the extension of the narrow-beam theory described in Section 2 to 
apply to a broad antenna system such as the SeaSonde, assuming ideal antenna patterns. From the 
antenna voltage cross spectra, we form as intermediate data products the first five Fourier angular 
coefficients of the broad-beam return over a selected range ring surrounding the radar. These 
coefficients, designated by the index n = -2, -1, 0, 1, 2, are defined in terms of the narrow-beam first 
and second-order return through the relation: 
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b);\cD) = \ a s 2 (co,<p)tf n (<p)d<p (29) 

where the integration over azimuth angle (p is performed over open water around the radar range 
cell and the superscripts refer to first- and second-order respectively. The narrow-beam radar cross 

1,2 

sections <J S (#>,#?) are defined in terms of the ocean wave spectrum by (15), (16). Following the 
notation in [5], the trigonometric functions $nW) a re given by 

tf n ((p) = sm(n(p)n<0 
= cos(n(p) n>0 

As described in [4], there are three steps in the interpretation of the radar spectrum to give deep- 
water wave information. 

a) The first- and second-order regions are separated. 

b) The first order region is analyzed to give the ocean wave spectrum at the Bragg wavenumber. It is 
assumed that deep-water theory is adequate for this step, as Bragg waves are short and hence 
insensitive to the effects f shallow water, see Figure 5. 

c) Second-order radar spectral data is collected from the four second-order sidebands of 10-minute 
averaged cross spectra and fit to a model of the deep-water ocean wave spectrum. Least-squares fitting 
to the radar Fourier coefficients is used to derive estimates of the significant wave height, centroid 
period and direction. During this step, the second-order spectrum is effectively normalized by the first- 
order, eliminating unknown multiplicative factors produced by antenna gains, path losses etc. 

Shallow- water analysis requires a further step: 

d) The shallow- water wave spectrum is calculated from the deep-water spectrum using (9)-(l 1). 

4.3 Model ocean wave spectrum 

For our analysis, we define a model for the deep-water ocean wave spectrum as the product of 
directional and nondirectional factors: 



S(k,q>) = Z(£)cos 4 \ J * j (31) 

The directional factor in (31) has a cardioid distribution around the dominant direction For 
describing the second-order spectrum, ^* is taken to be the dominant long-wave direction. For 
describing the first-order spectrum, ^ * is the short-wave direction, which is assumed to be the same as 
the wind direction. For the nondirectional spectrum we use the Pierson-Moskowitz model Z(k): 

a -0J4(k c /k) 2 

Z(k) = (32) 

whose parameters are the cutoff wavenumber ^ and a multiplicative constant A. The waveheight, 
centroid period and direction can be defined in terms of the model parameters. The significant 
waveheight follows from the directional spectrum through the relation: 



W = 4 



( y /2 

oo y 

j j S(k, (p)dkd(p 



(33) 
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This model has proven satisfactory for use in deep-water wave extraction software that produces 
waveheight, period, and direction. It has been used for real-time SeaSonde systems for many years, 
providing good agreement with in-situ measurements e.g. as shown in [4]. 

4.4 Results 

Figure 14 shows SeaSonde results in the second radar range cell calculated using the shallow-water 
theory described above, together with the ADCP results. 



Figure 14. SeaSonde (red) and ADCP (blue) results for (a) Significant waveheight (b) 
Wave period (c) Wave direction (d) Wind direction. 



■ 1 

o»: 1 












■ 
■ 
■ 
■ 
■ 
■ 
■ 


■ 
■ 

■ ■ 

■ ■ 

■ i 

— i—i — 

■ ■ 

■ ■ 

■ i 

■ ■ 

■ i 

■ i 
■ 






! 1 

r — * 
i 5" 


N 

* 




■ 
■ 

■ 
■ 

■ 
■ 
■ 

■ 
■ 


■ 

■ i 

■ i 

■ ■ 

■ ■ 

■ ■ 

: ,i t — J 




i i 
: ** 


i* i 

MM 


a— y — ^ 

x* 

"X ^ 




■ 
■ 
■ 
■ 
■ 


: iiS» rs sa 








■si 




■ ■ 

■ ■ 

■ ■ 

■ i 

■ i 

■ ■ 


■ ■ 

i ■ ■ 
i ■ ■ 
i ■ ■ 
i ■ ■ 
i ■ ■ 

■ ■ 






■ 
■ 
■ 



q 

mC 

> 

(T3 



02 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.C 

Days from start of December 29, 2005 




Sensors 2008, 8 



4632 



350 



300 



200 



160 



100 












































_ . 






1 



• 
■ 
■ 
■ 


i 
i 
i 


i 

i i 
i i 
i i 


?f *** 
--- 


flu* 
--- 




1* Ji 


jflMlfcil 




■ 
■ 
• 
■ 
■ 
■ 


■ 
■ 
■ 
■ 
■ 
■ 














■ 
■ 
■ 
■ 
■ 


■ 

■ 
■ 

j 













0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

Days from start of December 29, 2005 



1.8 2.0 



300 

c 
o 



200 



160 



100 



50 



1 






■ 


CJKKX" 






■ 


■ 
■ 
■ 
■ 
■ 


















— t— - 

■ 
■ 
■ 
■ 
■ 


















safe 




'A 


■ — 


■ 




















































0.2 0.4 0.0 0.8 1.0 1.2 1.4 1.0 

Days from start of December 29, 2005 



1.B 2.0 



It can be seen from Figure 14 that southerly winds veer to the northwest after the passage of a storm 
front. Subsequently the wave height and period increase suddenly. Spectral saturation may be 
occurring at the peak of the storm, causing overestimates in the waveheight. Wave direction remains 
about the same, as due to wave refraction, wave directions in very shallow water are nearly 
perpendicular to the depth contours. Both radar and ADCP are observing directions in shallow water, 
hence this perpendicular condition is being enforced on the longer waves, although the wind direction 
driving short waves is seen to change significantly over this storm period. 

Table 1 gives the bias and standard deviation between the SeaSonde and ADCP measurements of 
waveheight, wave period and direction, for the short-period waves before the storm and the longer- 
period waves afterwards. 
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Table 1. Comparison statistics, radar vs. ADCP. 



Before storm After storm 



Waveheight 

Standard deviation 0.25m 0.35m 

Bias -0.23m 0.17m 

Wave Period 

Standard deviation QJ6s 0.60s 

Bias 0.70s -0.27s 

Wave Direction 

Standard deviation 13.8° 19.7° 

Bias -9.5° 17.0° 



To emphasize the necessity of taking shallow water into account for this location, we estimated the 
waveheight assuming infinitely deep water. Figure 15 shows the SeaSonde results together with the 
ADCP waveheight. Clearly waveheight is overestimated with this assumption. The simulations 
described in Section 3 indicate that the cause of this overestimate is the failure to account for the 
increase of the coupling coefficient in shallow water. 

Figure 15. Significant waveheight: 

Red: SeaSonde calculated assuming infinite water depth. Blue: ADCP. 
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Table 2 gives the bias and standard deviation between the SeaSonde and ADCP waveheight 
measurements, with the former calculated assuming infinitely deep water. 

Table 2. Comparison statistics, radar vs. ADCP assuming deep water. 



Before storm After storm 



Waveheight 

Standard deviation 0.25m 0.95m 

Bias 0.19m 0.90m 



5. Conclusion 

We have presented the theory of narrow-beam HF radar sea-echo from shallow water and illustrated 
the effect of decreasing water depth using simulations for a simple swell model of the ocean wave 
spectrum. The second-order spectral energy increases relative to the first-order as the water depth 
decreases, resulting in spectral saturation when the waveheight exceeds a limit defined by the radar 
transmit frequency. This effect is particularly marked for lower radar transmit frequencies. For 
waveheights above the saturation limit, the perturbation expansions on which Barrick's equations (15), 
(16) are based fail to converge. The saturation limit on waveheight is less for shallow water. Shallow 
water affects second-order spectra (which gives wave information) far more than first-order (which 
gives information on current velocities). Figure 11 shows the depths at which shallow- water effects 
become significant plotted as a function of radar frequency for an lis wave. We discuss how the 
waveheight and depth limits would change for a more general model. 

The shallow-water theory was then extended to apply to broad-beam systems such as the SeaSonde 
and applied to the interpretation of two days of radar data measured by a 25Mhz SeaSonde located on 
the New Jersey shore. During the measurement period, a storm passed over the area. An ADCP was 
operated in the second radar range cell in water 8m deep. Radar results were compared with 
simultaneous ADCP measurements. The comparison confirms aspects of the theory presented in 
Section 3. For the longer period waves occurring after the passage of the storm front, the standard 
deviation between SeaSonde and ADCP waveheight measurements decreased by a factor of three when 
the effects of shallow water were included in the analysis, and the bias decreased by a factor of five. 
Possible explanations for the remaining discrepancies are (a) the assumption of parallel depth contours 
(b) the assumption that the wave spectrum is homogeneous in the circular radar range cell (c) 
saturation in the radar spectrum around the peak of the storm, which, as discussed in Section 3, leads to 
the over-prediction of the waveheight. 
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